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Abstract 

m _ 

. The first results from the KamLAND experiment have provided confirmational evidence 

I for the Large Mixing Angle (LMA) Mikheyev-Smirnov-Wolfenstein (MSW) solution to the 

solar neutrino problem. We do a global analysis of solar and the recently announced Kam- 
LAND data (both rate and spectrum) and investigate its effect on the allowed region in the 
I Am^ — tan^ 6 plane. The best-fit from a combined analysis which uses the KamLAND rate 

C3 ' plus global solar data comes at Am^ = 6.06 x 10~^ eV ^ and tan^ 6 = 0.42, very close to the 

global solar best-fit, leaving a large allowed region within the global solar LMA contour. The 
Qh, inclusion of the KamLAND spectral data in the global fit gives a best-fit Am2 = 7.17 x 10"^ 

{il^l eV and tan^ 6 = 0.43 and constrains the allowed areas within LMA, leaving essentially two 

^ • allowed zones. Maximal mixing though allowed by the KamLAND data alone is disfavored 

by the global solar data and remains disallowed at about 3a. The low Am^ solution (LOW) 
is now ruled out at about 5a with respect to the LMA solution. 
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1 Introduction 

It is fair to say that the recently announced first results of the Kamioka Liquid scintillator Anti- 
Neutrino Detector (KamLAND) experiment [1] constitute a highly anticipated milestone in our 
understanding and resolution of the three decade old solar neutrino problem. The origins of this 
puzzle lie in the early deficit measurements of the solar neutrino fiux in the pioneering Homestake 
chlorine experiment [2]. This discrepancy between the expected rate, as predicted by increasingly 
refined solar model calculations [3] and the measured one has been subsequently confirmed and 
buttressed over the years by results from the ''^Ga experiments SAGE, GALLEX and GNO [4, 5], 
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the Kamiokande and the Super-Kamiokande experiments (SK) [6], and most recently from the 
Sudbury Neutrino Observatory (SNO) [7, 8]. In particular, SK has provided valuable zenith angle 
and energy spectrum information in addition to total rate measurements of the high energy Boron 
flux, and SNO has provided crucial neutral current (NC) and charged current (CC) rate data 
along with spectrum results. Over the years, these experimental results have been culled together 
with our understanding of neutrino mass, mixing and resonant matter oscillations to obtain the 
allowed parameter space in terms of the mixing angle tan^ d and mass-squared difference Am^ of 
the neutrino states. The analysis of global solar data carried out by various groups favours the 
LMA solution based on MSW resonant matter oscillations [9] as the most probable resolution of 
the solar neutrino problem [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20]. 

KamLAND [21] is a 1 kton liquid scintillator neutrino detector, designed specifically to test 
the LMA solution. It is located at the earlier Kamiokande site in the Kamioka mine in Japan. 
Its main objective is to look for oscillation of i7f, coming from Japanese nuclear power reactors 
situated at distances ranging from ~ 80 km to 800 km. The bulk (~ 79%) of the measured flux 
is however from reactors which are at distances between 138 km to 214 km. The z/gS are detected 
via the inverse beta decay reaction t/g +P + n. Both the scintillation emitted by the positron 
as it moves through the detection medium, and its subsequent annihilation with an electron are 
recorded. The delayed coincidence of the positron with the 2.2 MeV 7-ray from the capture of the 
neutron constitutes a largely background free signal. The total visible energy corresponds 
to Ee+ +fne, where Ee+ is the total energy of the positron and rrie the electron mass. The positron 
energy is related to the incoming antineutrino energy as Ee+ = — E^ec — (^n — ^p) MeV 
{rrin — rrip = 1.293 MeV is the neutron-proton mass difference). Erec is the average neutron recoil 
energy calculated here using [22]. The energy resolution is a{E)/E = 7.5%/v^, E is in MeV. 

The first data from KamLAND gives the ratio of the observed number of events to the expected 
number of events to be [1] 

Rkl = 0.611 ± 0.085{stat) ± OMl{syst) (1) 

for an exposure of 162 ton-yr and a visible energy above 2.6 MeV ^. They have also presented the 
observed positron energy spectrum. 

In a pre-KamLAND analysis [23] we have shown that an energy integrated rate in the range 0.3 
-0.8 will provide confirmation for the LMA solution. In particular the solar LMA best-fit predicted 
a KamLAND rate of 0.65 which is close to the observed rate. This is the first confirmation of the 
LMA solution to the solar neutrino problem using terrestrial neutrino sources. We also showed 
that for a rate below 0.9 the LOW solution to the solar neutrino problem is disallowed at more 
than 3(7. Hence in this paper we focus on the LMA solution and perform a global analysis which 
combines 

(i) KamLAND rate and global solar data 

(ii) KamLAND spectrum and global solar data 

We find the allowed area from each of the above analyses and discuss their contributions in 
sharpening our knowledge of neutrino mass and mixing parameters. The current KamLAND 
and global solar data spht the allowed LMA region in two parts - a low Am^ region (low-LMA) 

^Below this energy the background due to the geophysical neutrinos dominate. 
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and a high Am? region (high-LMA), which has less (by ~ 2a) statistical significance. A more 
precise determination of Am^ and tan^ 6 should be possible with increased statistics and reduced 
systematics of the spectral data from KamLAND [24, 25, 26, 27, 28, 29, 30, 23]. We demonstrate 
the potential of 1 kton-yr spectral data in discriminating between the two allowed regions and 
further constraining the parameter values by simulating the spectrum at different values of Am^ 
and tan^ 9 selected from the allowed area of the global solar+KamLAND analysis. We find that 
if the true spectrum corresponds to that simulated at points in the low-LMA region then with 

1 kton-yr exposure the high-LMA part can be further disfavoured. For spectrum simulated at 
high-LMA values however the ambiguity between the two zones persists. 

2 Analysis and Results 

The total event-rate/sec in the KamLAND detector is given as [23] 

Nkl = j dE,a{E^)N,Y.S^?^^^^^ (2) 

where a{E^) denotes the cross-section; Si denotes the spectrum from a given reactor i and involves 
the neutrino spectrum from the fission of a particular isotope, the characteristic energy released 
per fission by the isotope and the fractional abundance of the isotopes. For further details of the 
spectrum, cross-section, fuel composition etc. we refer the reader to [23]. Np denotes the number of 
target protons. The declared KamLAND data corresponds to a fiducial mass of 408 ton, resulting 
in 3.46 X 10^^ free target protons[l]. Relative fission yields for the various fuel isotopes are also 
taken in accordance with [1], as is the integrated thermal power fiux of 254 Joule/cm^. Pi(Pe i^e) 
is the two-generation survival probability for the antincutrinos from each of the reactors i and di is 
the distance of reactor i to KamLAND in km. In addition, we include the event selection criteria 
used by the KamLAND collaboration corresponding to an efficiency of 78.3% [1]. 

We first do a statistical analysis of the KamLAND rate and global solar data. For KamLAND 
rate we define the ^ 

f r>expt r>theory\2 

Xkl — {■^) 

where a — '^JcTgyst + CTstat^ (^syst and astat being the total systematic and statistical error in the 
KamLAND data respectively (cf. Eq.(l)). 

Rkl = ^ (4) 

^^KL 

Nkl obtained from Eq.(2) with P(t'e I'e) = 1- 

For the solar analysis we define the function in the "covariance" approach as 

N 
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Figure 1: The 90%, 95%, 99% and 99.73%(3a) C.L. contours from a analysis using KamLAND 
rate + global solar data. The dashed line shows the presently allowed 3cr solar contour. 

where Ri are the solar data points, is the number of data points (80 in our case) and {(yfj)~^ 
is the inverse of the covariance matrix, containing the squares of the correlated and uncorrelated 
experimental and theoretical errors. We use the data on total rate from the CI experiment, the 
combined rate from the Ga experiments (SAGE+GALLEX+GNO), the 1496 day data on the SK 
zenith angle energy spectrum and the combined SNO day-night spectrum. For further details of 
our solar analysis we refer the reader to [12, 13]. 

The for the combined solar and KamLAND rate analysis is defined as 

X^ = xl + Xkl (6) 

The best-fit after including the KamLAND data comes at Am^ = 6.06 x 10"^ eV^ and tan^ 6 = 
0.42. Thus the best-fit point does not change significantly with respect to that obtained from only 
solar analysis [12]. In Figure 1 we draw the 90%, 95%, 99% and 99.73% C.L. allowed area in the 
LMA region from a combined solar+KamLAND rate analysis. Superimposed on that we show the 
3(T (99.73% C.L.) allowed area from solar data alone. Large area within the LMA regions is seen 
to remain allowed. 

Apart from the data on energy integrated total rates, KamLAND collaboration has also pro- 
vided the observed positron visible energy spectrum, albeit with low statistics. We incorporate 
this in our analysis to extract the shape information from this data as the spectral distortion is 
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Figure 2: The 90%, 95%, 
LAND spectrum data. 



99% and 99.73%(3a) C.L. contours from a analysis using the Kam- 



a very sensitive probe of Am? . For the KamLAND spectral data, we perform our analysis using 
a definition of x\ispec assuming the data to be Poisson-distributed which is appropriate for data 



with low statistics. For this case 



E 



Qexpt 

OfV Qiheory qexpt \ . c) Qexpt i / ^KL,i ^ 
^{J^nOxL,i ~ ^KL,i) + ■^^KL,i „ Mheory ' 



+ 



a 



(7) 



sys 



where asys is taken to be 6.42% [1] and a normalisation factor is allowed to vary freely. The sum 
is over the 13 KamLAND spectral bins. In [1] the errors for the shape distortion are attributed 
to energy scale, energy resolution, v^, spectrum and fiducial volume. A more refined statistical 
analysis would involve evaluating the systematic errors in each bin due to these sources at each 
Air? and tan^ 9 as well as taking into account of the background events and their errors in each 
bin. This will be possible as and when more detailed information will be available. 

For the spectrum analysis we get the best-fit values of Am? and tan^ 9 to be 7.17x 10~^ eV^ and 
0.64 respectively. This is close to that obtained by the KamLAND collaboration but our best-fit 
9 is not maximal as in [1] . Apart from the above there are other minimas with reduced statistical 
significance. In Table 1 we present the best-fit values and Xmm t^e global minima and the 
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Figure 3: The 90%, 95%, 99% and 99.73% C.L. contours from a analysis using KamLAND 
spectrum data along with global solar data. The dashed line shows the presently allowed 3a only 
solar contour. 

second minima which is obtained at a higher Am^ value. In Figure 2 we present the allowed 
areas in Am^ — tan^ 6 plane from KamLAND spectrum analysis. This is understandably slightly 
different from what KamLAND has obtained in [1]. The KamLAND data analysis procedure as 
outlined in [1] is somewhat different and the full details are not known to us. 

Next we do a combined analysis of KamLAND spectral data together with the global solar 
data. The the combined analysis is defined as the sum of the individual contributions. We 

present the results in Table 1. The best-fit comes at Am^ = 7.17 x 10~^ eV^ and tan^ 6' = 0.43. 
From Table 1 we also see that there is a second minima at a higher Am^ with a reduced statistical 
significance by about 2a with respect to (w.r.t.) the global minima. 

In Figure 3 we show the combined allowed area with the solar and KamLAND spectrum 
analysis superimposed on the 3a solar contour. The allowed area is seen to be much constricted 
with the inclusion of the KamLAND spectral data. At 90% C.L. only a small region about the 
best-fit point remains allowed. At 99% C.L. however there are two distinct allowed zones - one 
around the global best-fit point (low-LMA) and the other around the higher Am^ corresponding to 
the second minima (high-LMA). The former is preferred by the KamLAND data and to a greater 
extent by the global solar data. At 99.73% the demarkation between the two zones disappear. 

In Table 2 we show the allowed ranges of the values of the parameters at 99% C.L. obtained 
from the global analysis including the solar and KamLAND spectrum data. The allowed range 
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KamLAND 


7.17 X 10- 


-5 


0.64 


5.71 


1.50 X 10- 


-4 


0.34 


8.24 


KamLAND + Solar 


7.17 X 10- 


-5 


0.43 


74.39 


1.48 X 10" 


-4 


0.44 


81.51 



Table 1: The Xmin best-fit values of the oscillation parameters obtained from the analysis 

of the KamLAND spectrum data alone and from the global analysis of the KamLAND spectrum 
data and the solar neutrino data. 



Allowed 
Zone 


99% C.L. Range of 
Am? in eV^ 


99% C.L. Range of 
tan^e 


low-LMA 
high-LMA 


5.3 X 10-5 < Am^ < 9.9 x 10"^ 
1.3 X 10-^ < Am^ < 1.8 X 10"^ 


0.28 < tan^ 9 < 0.79 
0.34 < tan^ < 0.55 



Table 2: Range of parameter values allowed at 99% C.L. from the global solar and KamLAND spec- 
trum analysis. 



of tan^ 6 is not reduced in the preferred low-LMA zone with the inclusion of KamLAND data, 
although in the high-LMA zone it is somewhat restricted. 

For the LOW solution we get Xmm — 97.28. This implies that LOW is now ruled out at 4.4(T 
w.r.t. the LMA solution. The maximal mixing solution {Xmin — 88.89) is disfavored at 3.4(7 (w.r.t. 
LMA) from the combined KamLAND -|- solar analysis. 

3 Projected Analysis 

In this section we explore whether the future KamLAND spectrum data would be able to de- 
termine the allowed zones more precisely. In particular, it would be interesting to see whether 
KamLAND , either by itself or in conjunction with the solar data, can choose between one of the 
two allowed islands, and pin down the values of the mass and mixing parameters unambiguously. 
We try to look into this by doing a statistical analysis of the projected 1 kton-yr KamLAND data. 
We choose sample values of Am^ and tan^ 6, from the allowed zones obtained from the global 
solar -l-KamLAND analysis (cf. fig. 3) and simulate the spectrum at these points for 1 kton-yr of 
data (approximately 2.5 years of KamLAND live-time ), using a randomizing procedure which 
takes care of the fluctuations. We use these simulated spectra in a analysis and reconstruct 
the allowed regions in the Am^-tan^ 9 plane. In Figure 4 we display these reconstructed allowed 
regions from 1 kton-yr projected KamLAND spectrum alone. Figure 5 shows the 1 kton-yr sim- 
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Figure 4: The 90%, 95%, 99% and 99.73% C.L. contours for using the 1 kton-yr projected Kam- 
LAND spectrum. The different panels are for the simulated spectrum at values of Am^ and tan^ 9 
indicated by the black dots. The dashed line shows the currently allowed 99% C.L. contour from 
solar+KamLAND analysis. 

ulated spectrum with errorbars at the Am^ and tan^ 9 corresponding to each of the panels of 
Figure 4. In Figure 6 we present the allowed regions from a combined analysis of global solar data 
and 1 kton-yr KamLAND spectrum data, simulated at the same set of points as in the previous 
two figures. The dashed lines in Figures 4 and 6 give the current 99% C.L. allowed regions. For 
the analysis with the projected KamLAND spectra we assume Gaussian statistics, which is more 
appropriate in this case. 

First let us make a note on the rationale of the representative values of Am? and tan^ 9 chosen to 
simulate the spectrum. Panel 1 corresponds to the KamLAND spectrum simulated at the low-LMA 
best-fit while panel 7 corresponds to that generated at the high-LMA best-fit. These are the two 
favored points from the current data. We also simulate the 1 kton-yr KamLAND spectrum at few 
other points deviated from the best-fits. In panel 2(3) we have chosen Am^ corresponding to the 
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Figure 5: The 1 kton-yr simulated KamLAND spectrum for the different sets of Am^ and tan^ 9 
corresponding to Figures 4 and 6. The histogram shows the unoscillated spectrum for 1 kton-yr. 

low-LMA best-fit but a lower(higher) value of 9. Similarly panel 8 is for Am^ corresponding to the 
high-LMA best fit but at a higher value of 9. The panels with same Am^ but different 9 are chosen 
to demonstrate the impact of 9 on the reconstructed regions. We choose Am? values lower(higher) 
than the low(high)-LMA best-fit in panels 4(5). Panels 6 and 8 display the reconstructed regions 
for Am^ values lower and higher than the high-LMA best-fit respectively. These set of values give 
adequate coverage for studying the projected sensitivity of the reconstructed regions on the choice 
of Am^ and 9 currently allowed at 99% C.L.. 

Figure 4 addresses the issue if future KamLAND data can by itself make a sharper demarkation 
between the two allowed islands and the dependence of this on the oscillation parameters. The 
figure shows that if we choose the simulation point at higher Am^ values (panels 5, 6, 7, 8 and 9) 
then there are large allowed regions from the spectrum data, extending upto 10^'^ eV^. The Figure 
5 shows that for these values of Am^ the spectral suppression tends to become flat (undistorted) 
leading to increased fuzziness. The same is also true for lower values of 9. The lesser distortion at 
low values of tan^ 9 is due to a diminished energy dependence resulting from a small oscillatory 
term in the survival probability [23]. This also leads to an increase in allowed area (cf. panel 
2) although to a lesser extent as compared to the panels with high Am^ values discussed earlier. 
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Figure 6: Same as in Figure. 4 but from a combined analysis using the solar+CHOOZ+1 kton-yr 
projected KamLAND spectrum data. Since the solar data disfavour the 6* > 7r/4 regions this 
figure is plotted with tan^ 6* extending upto 1.0. 

There are also some allowed regions at Am^ values lower than that allowed by the current 99% 
CL. contour in panels 5, 6, 7 and 9. However the figure reveals that the allowed Am^ range 
around the low-LMA and high-LMA zones get reduced in most of the panels. With 1 kton-yr 
KamLAND spectrum data the region around Am? ~ 10^'* eV^ present in Figure 3 gets disallowed 
in all the panels and the low-LMA and high-LMA regions get bifurcated even at 3cr level. In 
general, tighter constraints in the Figure 4 are associated with spectra with a shape significantly 
different from the no oscillation spectrum. For the KamLAND baselines this corresponds to a 
lower Am? than the present best-fit. Panel 4, which is at the global solar best-fit, shows that 
for such cases just the 1 kton-yr spectrum data from KamLAND can pick out a sharply defined 
allowed zone around the simulation point unambiguously. 

In Figure 6 we show the CL. allowed regions from a combined analysis of the global solar + 
CHOOZ + 1 kton-yr KamLAND simulated spectrum data. Through these plots we investigate the 
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impact of the current solar data to resolve the ambiguity still admitted by the 1 kton-yr projected 
KamLAND spectrum data. We find that the solar data is instrumental in ruling out a large part 
of the parameter space allowed by the 1 kton-yr KamLAND only analysis. A comparison of the 
panels 1, 2, 3 and 5 in the figures 4 and 6 shows that for spectrum simulated in the low-LMA 
region the inclusion of the solar data reduces the statistical significance of the high-LMA zone. It 
also disallows the high Am^ regions above 2 x 10~^ eV^. In these regions the solar data requires 
a fiux normahsation factor (/s) ~ 0.8 which is in confiict with the SNO NC measurement of 
/b ~ 1.0, thus disfavoring these zones. For spectrum simulated in the high-LMA region in panels 
6 to 9, even though the large allowed regions beyond 2 x 10^"^ eV^ as well as the allowed regions at 
low Am^ get mostly removed by the solar data, the ambiguity between the allowed islands remain. 
In fact a comparison with Figure 4 reveals that for these panels the inclusion of the solar data 
increases the statistical significance of the low-LMA allowed regions. The KamLAND spectral data 
allows the low-LMA region only at 99% C.L. in panels 7, 8 and 9, and at 99.73% C.L. in panel 6. 
However with the inclusion of the solar data in the analysis, these regions become allowed at 90% 
and 99% C.L. respectively, as the solar data prefers the low-LMA zone. Therefore to remove the 
ambiguity for the spectrum corresponding to the high-LMA zones, one would require KamLAND 
data with higher statistics, which will be able to determine the spectral shape and hence Am^ 
more precisely. 

4 Conclusions 

In conclusion, we have investigated the impact of the first results from KamLAND on neutrino 
mass and mixing parameters in conjunction with the global solar neutrino data. KamLAND is 
completely consistent with the LMA solution, to the extent that the observed KamLAND rate is 
close to that predicted by the best-fit point of the LMA solution to the solar neutrino problem. 
As a result the combined analysis of the solar and KamLAND rates data allows a large area 
within the solar LMA region and the global solar best-fit does not change much with inclusion of 
KamLAND rates. The constraining capabilities of the spectrum data is much stronger and with 
only 145 days observed spectrum KamLAND can exclude certain parts of the LMA parameter 
space. After including the spectral data the allowed LMA zone consists mainly of two disconnected 
regions, one around the best-fit and another at a higher Am^. The two zones merge at 3(t. Maximal 
mixing though allowed by the KamLAND alone, is found to be still disfavored by the combined 
solar and KamLAND data at more than 3o". The LOW solution which was allowed at 3o" from the 
global solar data and which predicts null oscillations in KamLAND is now disfavored at almost 
5(7 w.r.t the LMA solution. 

With LMA now confirmed, the next focus of KamLAND would be a more accurate deter- 
mination of the mass parameter by distinguishing between the two allowed sectors in the LMA 
region. We have explored this through a projected analysis with 1 ktyr simulated data. With 1 
kton-yr projected spectrum data the allowed Am^ ranges around both low-LMA and high-LMA 
zones decrease in size and they get separated at 3a by the spectrum data itself. The inclusion 
of the solar data disfavours (favours) the high(low)-LMA zone if the spectrum is simulated in the 
low(high)-LMA area. Thus the allowed areas become more precise for low-LMA spectrum while 
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ambiguity between the two zones remains for high-LMA spectrum. A higher statistics from Kam- 
LAND is expected to resolve this ambiguity. A more precise determination of the mixing angle 
will be possible from a more accurate measurement of the CC/NC ratio at SNO. 
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